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In the context of bipartite bosonic systems, two notions of classicality of correlations can be 
defined: P-classicality, based on the properties of the Glauber-Sudarshan P-function; and C- 
classicality, based on the entropic quantum discord. It has been shown that these two notions 
are maximally inequivalent in a static (metric) sense — as they coincide only on a set of states 
of zero measure. We extend and reinforce quantitatively this inequivalence by addressing the dy¬ 
namical relation between these types of non-classicality in a paradigmatic quantum-optical setting: 
the linear mixing at a beam splitter of a single-mode Gaussian state with a thermal reference 
state. Specifically, we show that almost all P-classical input states generate outputs that are not 
C-classical. Indeed, for the case of zero thermal reference photons, the more P-classical resources 
at the input the less C-classicality at the output. In addition, we show that the P-classicality at 
the input — as quantified by the non-classical depth — does instead determine quantitatively the 
potential of generating output entanglement. This endows the non-classical depth with a new op¬ 
erational interpretation: it gives the maximum number of thermal reference photons that can be 
mixed at a beam splitter without destroying the output entanglement. 

PACS numbers: 03.67.Mn, 42.50.Dv 


I. INTRODUCTION 


Since the early days of quantum mechanics consider¬ 
able efforts have been spent in establishing whether a 
given physical system possesses genuinely quantum fea¬ 
tures. As far as bosonic systems are concerned, Wigner 
first attacked this problem introducing a quantum ana¬ 
logue to the classical phase-space [1]. Later on, a system¬ 
atic approach was finally developed in the framework of 
quantum optics, with the introduction of various classes 
of quasi-probability distributions defined over the quan¬ 
tum phase-space. Specifically, the analytical features 
of such distributions unveil physical constraints: when- 
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ever the normally-ordered distribution function — called 
Glauber-Sudarshan P-function [2, 3]— behaves like a reg¬ 
ular probability distribution, the corresponding state can 
be described as a statistical ensemble of classical fields 
and, in this sense, it cannot show any non-classical fea¬ 
ture [4]. In the following, these states will be referred to 
as P-classical states. 

On the other hand, the more recent development of 
quantum information theory promoted a reconsidera¬ 
tion of the quantumness of physical systems from an 
information-theoretical perspective. Since quantum sys¬ 
tems can be correlated in ways unaccessible to classical 
ones, the discrimination between classical and nonclas- 
sical states of a given system is pursued by studying 
the nature of the correlations among its subparts. In 
particular, quantum entanglement accounts for quantum 
correlations that may lead to the violation of local real¬ 
ism [5]. Moreover, even separable (i.e., non-entangled) 
states have been recognized to retain non-classical fea¬ 
tures, leading to the introduction of an entropic measure, 
called quantum discord, to capture the quantum features 
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of correlations beyond entanglement [6]. Following this 
criterion, classical states can be defined as states with 
vanishing discord, and we will refer to this notion as C- 
classicality. 

Although both acceptable and well-grounded, the two 
foregoing notions of classicality have been shown to be 
radically different, indeed maximally inequivalent, in the 
following sense: only a zero-measure set of states is clas¬ 
sical according to both criteria [7]. Besides embodying 
a matter of fundamental interest, this conclusion is also 
relevant for practical purposes, since it enlightens dif¬ 
ferent resources in quantum information processing [8]. 
However, such a characterization is based on purely ge¬ 
ometrical considerations and, as a consequence, it is in¬ 
trinsically static. In particular, the relation between C- 
and P-classicality in common physical processes remains 
unclear. In addition, a quantitative comparison between 
these two notions in terms of their respective figures of 
merit still lacks. 

In this work, we address the above issues in the con¬ 
text of quantum optics, whose description in term of 
the phase-space offers a natural framework to develop 
a quantitative analysis [9]. A paradigmatic setting in 
quantum optics is constituted by Gaussian states and op¬ 
erations, due to their relevance for quantum technologies 
and their thorough theoretical characterization [10-12]. 
Specifically, we address the dynamical relation of P- and 
C-nonclassical states arising from the linear mixing of 
Gaussian states at a beam splitter. In this setting, the 
availability of analytical expressions to quantify Gaus¬ 
sian P-classicality — in terms of the non-classical depth 
— and Gaussian discord and entanglement, is crucial to 
work out a quantitative comparison between the various 
notions of non-classicality. 

In particular, we consider the mixing of a generic Gaus¬ 
sian state with a reference thermal state, and explore the 
relationships between the classicality of the input state 
and the P- and C-classicality at the output. Specifically, 
we show that almost all P-classical input states gener¬ 
ate output states that are not C-classical. Indeed, for 
the case of zero thermal reference photons, the more P- 
classical resources at the input the less C-classicality at 
the output. These findings strengthen the inequivalence 
between P- and C-classicality by quantitatively extend¬ 
ing it to a process in which correlations are dynamically 
generated, rather then statically analyzed as in Ref. [7]. 
In addition, we show that the P-classicality at the in¬ 
put does instead determine quantitatively the potential 
of generating output entanglement. This endows the non- 
classical depth with a new operational interpretation: it 
gives the maximum number of thermal reference photons 
that can be mixed at a beam splitter without destroying 
the output entanglement. 

The paper is structured as follows. In Section II we 
give a brief account on Gaussian states and their phase- 
space representation, focusing on their bilinear interac¬ 
tion in linear optical devices. In Section III we review the 
two notions of nonclassicality and establish the notation 


for Gaussian discord, non-classical depth and entangle¬ 
ment used in the following. The reader familiar with the 
foregoing topic can skip the respective sections. In Sec¬ 
tion IV we analyze in details the generation of P- and 
G-nonclassicality by mixing of a Gaussian state with the 
vacuum, whereas in Section V we focus attention to the 
mixing with a thermal state, also introducing the concept 
of effective nonclassicality. Section VI closes the paper 
with some concluding remarks. 


II. LINEAR MIXING OF GAUSSIAN STATES 

The simplest bilinear interaction involving two bosonic 
field modes described by the annihilation operators di 
and 0,2 (with [ofe, 6,^] = 1) corresponds to the mode mix¬ 
ing and it is described by an effective Hamiltonian of 
the form Hj (x (a|a 2 -f oidl). This kind of interaction 
is very common in different quantum systems, ranging 
from optical modes in linear optical devices [13] to col¬ 
lective modes in ultracold atoms [14], opto- and nano¬ 
mechanical oscillators [15-18] and superconducting res¬ 
onators [19, 20]. For the sake of clarity, in this paper we 
focus on the quantum optics realm and we address the 
correlations properties of the two optical modes emerg¬ 
ing from a beam splitter (BS) when the input ones are 
excited in Gaussian states. 

Gaussian states (GSs) are states with Gaussian Wigner 
functions [21] and an exhaustive information about them 
is provided by the knowledge of the first and second sta¬ 
tistical moments of the quadrature operators. Informa¬ 
tion about correlations is contained in the second mo¬ 
ment and from now on we set the first moments to zero, 
without loss of generality. Upon introducing the quadra¬ 
ture operators g = (a + a^)/\/2 and p = {d — d^)/{i^/2), 
the covariance matrix (CM) of a single-mode GS of 
d is defined as [cr]^, = ^{{Rk,Ri}) - {Rk){Ri), with 
k,l = 1,2, being = (Pi,P 2 ) = the vec¬ 

tor of the quadratures and {•,•} the anticommutator. 
Now, the canonical commutation relations take the form 
[Rk, Ri] = iuJki, where ujki = (l~<5fe/)(~l)* are the entries 
of the is the 2x2 symplectic form u>. The set of the eigen¬ 
values {q,p) € of the position and momentum-like op¬ 
erators, endowed with the symplectic form a;, spans the 
real symplectic space F = (M^, u>), which is referred to as 
the phase space of the mode di. A single-mode GS may 
always be written as 

Qi{nt,ns) = S{r)iy{nt)S\r), (1) 

where S{r) = exp [^ (ra'^^ — r*a^)], r € C, is the squeez¬ 
ing operator and 

n{nt) = {nt + 1)“^ [nt/{nt + 1)]“^“ , 

is a thermal state with nt average number of photons, the 
quantity Ug = sinh^ jrj will be referred to as the number 
of squeezed photons. 
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By choosing a suitable rotating frame, the lossless BS 
exchange interaction is described by the unitary evolu¬ 
tion Ut = exp[d(d|d 2 — dia\)], with 0 € K and where 
T = cos^ 0 denotes the transmissivity of the BS. If 
T = 1/2, the BS is said to be balanced. Being a bilinear 
interaction of modes, this evolution preserves the Gaus¬ 
sian character of the state, and in turn induces a sym- 
plectic transformation Sr in the quantum phase space of 
the composite system, namely 


Sr 


\/t1 Vl — r I \ 

— Vl — T I ^/Tl J 


( 2 ) 


where I = diag(l,l). Given two uncorrelated single¬ 
mode GSs, with GMs cri and (T 2 , respectively, the sym- 
plectic transformation Sr, acting by congruence on the 
initial CM So = <ti 0 (T2, leads to the evolved CM 
S = Sr So Sr. As mentioned above and schematically 
depicted in Fig. 1, in this work we will consider a bi¬ 
partite quantum system of modes Si and 0 - 2 - The mode 
Si is initially in the zero-mean GS g = g{ns,nt) while 
mode S 2 is in a thermal state ly = v{n 2 )- Without loss 
of generality we will assume a real squeezing parameter 
r S K. Using this parametrization, the average number 
of photons in the hrst mode (S|Si)g = Tr[S|Si g{ns,nt)] 
explicitly reads: 

{d\di)g = Tit + {1 + 2nt)ns = ni- (3) 



FIG. 1: (Color online) Linear mixing of Gaussian states. The 
two input modes Si and S 2 , initially excited in the zero-mean 
Gaussian state g = g{ns,nt) and in the thermal state v = 
1 /( 112 ) respectively, enter a beam splitter of transmissivity r, 
after which quantum correlations are eventually established. 


In the phase space, the GSs g and 1 / are represented 
by the 2x2 GMs: 


1 


(Tg = diag - + ni + A, - 0 ni - A 


(4) 


and 


= [-+712 ]I, 


(5) 


respectively, where A = (1 + 2nt)\/nsi^ + ris)- Since 
the initial state Rq of the bipartite system is chosen to 
be factorized, namely Rq = g ® v, the total number of 
excitations is given by (ajoi + d 2 fl 2 )ito = ni + n 2 = N. 
Now we let the state i?o evolve through a lossless BS of 
transmissivity t using the symplectic transformation 
given by Eq. (2). The 4x4 CM of the two-mode output 
state R = UrRoUl is S = Sr{(Tg 0 cr,y)S^ and it reads: 


S = 


a+ 0 c+ 0 
0 a_ 0 c_ 
c+ 0 64 . 0 

0 c_ 0 5_ 


where: 
a± = 

b± = 

c± = 


^ + «2 j (1 - r) 0 Q 0 ni ± A ) T , 


^ 0 712 j T + Q + ni ± A ) (1 - r), 


( 6 ) 


1 


712 


1 


0 Til ± A 


(7a) 
(7b) 
(7c) 


It is useful to introduce the following local symplectic in¬ 
variants, which in this case are given by Ii = a_|_a_, I 2 = 
b^b -,13 = c+c_ and I 4 = detS. Via symplectic di- 
agonalization (tab can be cast into the diagonal form 
diag(A+, A+, A_, A_) where the expression of the sym¬ 
plectic eigenvalues is given by [ 22 ]: 


//i 0 /2 0 2/3 ± v/prT 72 0 2 / 3)2 


Positivity of gAB requires A_ > 1/2. 


III. NONCLASSICALITY FOR BOSONIC 
SYSTEMS 

Let us now review in some details the concepts of non- 
classicality we are going to consider, together with their 
respective figures of merit. The reader familiar with those 
concepts can skip this Section. 


A. Nonclassicality in the phase space: P-classicality 

Any bipartite bosonic state described by the density 
matrix gAB can be always expanded in terms of coherent 
states as follows: 

QAB = [ (fa[ d^/3P(a,/3) |a) (a| 0 1^) (/3| , (9) 

Jc Jc 

where la) and |/3) are coherent states of the two modes 
and P(a,/3) is the Glauber-Sudarshan P-representation 
of the state. P(a,/3) provides a complete characteriza¬ 
tion of the state. Equation (9) suggests that the state 
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of the electromagnetic field can be regarded as a mix¬ 
ture of coherent states, weighted by P{a,fi). However, 
in general the P-function cannot be regarded as a proba¬ 
bility density function. On the other hand, when all the 
conditions for the P-function to be a probability density 
are satisfied, one can conclude that it describes a clas¬ 
sical state of the bosonic field, motivating the following 
definition: 


Definition 1 (P-classicality) A state qab of a two-mode 
bosonic field is called P-classical if P{a,j3) is a regular 
and normalized positive function. 


In the case of a single-mode state g, we can introduce 
a generalized s-ordered Wigner function which encom¬ 
passes all the quasi-probability distributions : 

Wsia)= [ ( 10 ) 

Jc 

where D{a) = exp(Q;a^ — a*a) is the displacement oper¬ 
ator. In the case of s = — 1, 0,1 one recovers the Husimi, 
Wigner, and Glauber-Sudarshan functions, respectively. 
The latter, more than any other representation, can de¬ 
part from being a well-behaved probability density. In 
order to understand this fact, let us observe that the 
s-ordered Wigner function of a state is related to the P- 
function (s = 1) of the same state through a Gaussian 
convolution, namely: 


Ws{a) 


2 

7r(l — s) 


d^/3 exp 


2|a-/3p 

1 — s 


PiP), 

( 11 ) 


that can be seen as a smoothing operation. Given the 
P-function of the state of interest, as the parameter s 
moves towards —1, the resulting distributions Ws{a) get 
smoother and smoother. Since for s = — 1 the Husimi 
Q-function is recovered, we obtain a continuous interpo¬ 
lation between P- and Q-function, and we are guaranteed 
that this smoothing operation always succeeds in giving 
a true probability distribution. 

Based on this, it is possible to define a quantitative 
measure of P-nonclassicality, that is the non-classical 
depth of a quantum state [23, 24]. To this aim it is 
useful to introduce the parameter T = (1 — s)/2. For 
T large enough, Eq. (11) leads to a P-classical state and 
the smoothing operation is referred to as complete. If 
denotes the set of all T which give a complete smoothing 
of the initial the P-function, the non-classical depth is 
defined as 


T„=myT). (12) 


The non-classical depth ranges from 0 for coherent states 
to 1 for Fock states [25], whereas for a single-mode GS g 
we have [26]: 


Tm = max 


1 

2 




(13) 


where p, = Tr[^p^] is the purity of the state and r is the 
squeezing parameter introduced in Eq. (1). 


B. Nonclassicality of correlations: C-classicality 

The total amount of correlations between two classi¬ 
cal systems A and B is quantified by the mutual in¬ 
formation I{A : B) = H{A) -\- H{B) — P[{A,B) where 
H(X) = —^^pxix)logpxix) is the Shannon entropy 
of the random variable X. By exploiting the rela¬ 
tion pAB(P',h) = PA{o\b)pB{b>) one also gets the equiv¬ 
alent expression of the mutual information I {A : B) = 
H{A) — H{A\B) in terms of the conditional entropy 
hIa\B) = -Y,a,bPB{b) PA{a\h)\ogpA{a\h). The first 
expression of the mutual information has an immedi¬ 
ate extension to quantum systems, simply by replacing 
the Shannon entropy with the Von Neumann entropy 
S[g] = -Tr[eloge], namely/ m(£'ab) = 5'(pa) + <S'(pb)- 
S{gAB), if we address GSs, it can be expressed as [22] 

iMigAB) = f {^/h) +f (^ 72 ) -/ (A+)-/ (A_), (14) 

where f{x) = {x-\- P) log(a; -\-\) — {x—\) log(a; — 1). On 
the other hand, the extension of the second expression to 
the quantum realm involves a measurement on one of the 
two parties, say P, described by the positive operation¬ 
valued measure (POVM) {Hfe}, H^ > 0, = P 

The probability to obtain the outcome k is given, ac¬ 
cording to the Born rule, by pk = Tr[£»^B I G H^,] and 
the conditional state of A with respect to the k outcome 
g^'^Q = (pfc)~^Ti'B[t'yiB I G Hfc]. The maximum amount 
of information we can gain on the part A by locally mea¬ 
suring the other part thus has the non-trivial expression 

Ca\b{0ab) = 1*^ (pa) - ^PkS I , (15) 

and invokes an optimization procedure over the set of all 
measurements. The quantum discord is properly defined 
as the difference between these two quantities [27] 

PAlBiSAB) = ImigAB) - Ca\b{0Ab) ■ (16) 

We can thus conclude that a system shows some quan¬ 
tumness as soon as the discord is different from zero, 
providing us with the following criterion: 

Definition 2 (C-classicality) A state gAB of a two¬ 
mode bosonic field is called C-classical if Paib^Pab) = 
'Pb\a{0ab) = 0. 

If we restrict to the subclass of GSs and Gaussian mea¬ 
surements, an analytical expression of the quantum dis¬ 
cord can be derived, which is called Gaussian discord, 
and is given by 

PMBiPAB) = f (/^) +/ (/f2) - / (A+) - / (A_) , 

(17) 

where P™)) an analytical expression as a function 
of the local symplectic invariants Ik [28-32]. It is worth 
noting that for a large class of Gaussian states [33] the 
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Gaussian discord of Eq. (17) coincides with the quantum 
discord, i.e. the maximum in Eq. (15) is achieved by a 
Gaussian measurement. 

A stronger form of quantum correlations with respect 
to discord is given by quantum entanglement. In the case 
of two-mode GSs a necessary and sufficient condition can 
be derived to assess the presence of entanglement [34]. It 
is essentially based on the positivity of qab under partial 
transposition (PPT), that is the positivity of the density 
matrix obtained by the transposition applied only to one 
part of a system [35]. One can show in fact that the 
symplectic eigenvalues of the partially transposed states 
are given by: 


It is apparent that the effect of the evolution only 
amounts to a re-parametrization of the argument, that 
does not affect the functional form of both Pg and 
Po- Thus we can conclude that output state R is P- 
nonclassical if and only if g is P-nonclassical, that is, for 
our configuration the two-mode P-nonclassicality of the 
output equals single-mode P—nonclassicality of the input. 
Let us remark, for the sake of Sec.V, that the foregoing 
argument applies as well for the case of mixing with a 
reference thermal state with positive temperature, given 
that the P-function of a general thermal state is a well- 
behaved probability density. The non-classical depth in 
Eq. (12) relative to the mode di, can be expressed as: 


A+ — 


'h+h- 2/3 ± ^ih +I2- 2/3)2 - 4/4 


, ( 18 ) 


Tm = max 


1 - 2 u 
2 



( 21 ) 


and a Gaussian state qab is entangled if and only if 
A_ < 1/2. In fact, a measure of entanglement is 
given by the logarithmic negativity [36], that is E{cr) = 


max 


log(2A_),0 


where u = 1 -I- ni — A is the minimum eigenvalue of 
the CM (Tg, as it is apparent from Eq. (4). The condi¬ 
tion Tm = 0 singles out a P-classicality threshold, which 
can be made explicit either as a function of the thermal 
component n*, hence having 


IV. NON-CLASSICALITY ARISING FROM 
MIXING A GAUSSIAN STATE WITH THE 
VACUUM 


n 


p 

S 


nt 


1 -|- 2nt 


or of the squeezed ones 


( 22 ) 


We start considering the case in which the reference 
input state of mode 0,2 is the vacuum, i.e. 712 = 0 and, in 
particular ctj^ —> (Tq = and ni = N. 

A. P-classicality 

The initial state of the system Rq is clearly C-classical. 
On the contrary, P-nonclassicality has to be addressed 
both in the input and output channels, in order to see 
wether differences arise. The Glauber-Sudarshan P- 
representation of Pq is given by 

Ro= f (fa [ fp Pg{a) PfP) |a) (a| 0 |/3) (/?| , (19) 
Jc Jc 

where the P-function Pr^ is factorized in the product 
of Pg and Po, being the latter the Glauber-Sudarshan 
P-function of the vacuum. Moreover, since Pq is a well- 
behaved probability density, any possible P-nonclassical 
feature of the input state is due to the pathological be¬ 
havior of Pg alone, and is quantified by the nonclassical 
depth of Eq. (12). Since the action of the BS evolution 
on the two-mode displacement operator is Ur Da{a) G 
Db{/3)Ul = Daiy/ra -\- ^/1 - TP)Db{y/ffi - Vl - Ta), 
and it amounts to a rotation of the arguments, one ob¬ 
tains the following P-representation of R: 

PR{a,f3) = Pgi^/ra - Vl - rjd) Po(Vl - ra -\- ^/rfi). 

( 20 ) 


n^^ = ns + \/ns{l + ns). (23) 

Whenever the average number of squeezed photons ex¬ 
ceeds or the thermal component falls below n/ the 
state Q turns out to be P-nonclassical. From now on, 
otherwise differently stated, shall be employed as the 
P-classicality threshold. 


B. C-classicality and generation of Gaussian 
discord 

We now address general quantum correlations, and in¬ 
vestigate the generation of Gaussian discord. Being the 
parties involved the output modes bi = W diU and 
62 = 02 U, we will refer to Pi |2 as the 61 -discord 

and P 2 I 1 to the 62 -discord; when they both coincide the 
symbol P will be employed employed. In fact, this is 
the case of a balanced BS, namely Pi| 2 (ns,nt, 1/2) = 

T’2|i(ns, nt, 1/2) \/Us, nt. 

In Fig. 3 we show a plot of the Gaussian discord as a 
function of the squeezed and thermal component of the 
input state g, for the balanced case r = 1/2. Except for 
the trivial case of a vacuum input state g = | 0 ) ( 0 |, it is 
apparent that the discord is always positive and there¬ 
fore, contrary to the case of P-classicality, there is no 
C-classicality threshold: whatever the input state, a bal¬ 
anced BS is capable of generating quantum correlations. 
This is also in agreement with the fact that typically al¬ 
most all states possess positive discord [45]. 
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FIG. 2: (Color online) Nonclassicality by mixing with the 
vacuum: Plot of the implicit function Tm = 0, which defines 
the P-classicality threshold. It takes the explicit expression 
of Eq. (22) as a function of nt, while as a function of Us it is 
given by Eq. (23). In this case the non-classicality threshold 
and the separability threshold nf’’ coincide. The black 
lines are curves of fixed total energy N = Us + nt + 2nsnt ■ 
dashed for odd values N = 2fc+l, A: = 0,1,... and dot-dashed 
for even values N = 2k, fc = 1, 2,... of the total energy. 


From Fig. 3 a quantitatively relevant feature emerges. 
Considering input states below the P-classicality thresh¬ 
old (denoted by the dashed black curve), one can observe 
that the output discord and hence the C-nonclassicality 
increases as the input non-elassical resources decreases. 
This is true regardless the constraints that one considers: 
either moving along the curves at constant Us, nt, or to¬ 
tal energy N, the discord increases as Hs decreases or n* 
increases. This is a quantitative feature that strikingly 
confirms — together with the absence of a C-classicality 
threshold — the inequivalence between the two notions 
of classicality considered here. 

In Fig. 3 the Gaussian discord corresponding to three 
families of input states has been highlighted: the discord 
generated by a thermal input state = V{Q,nt) cor¬ 
responds to the solid blue curve; the dotted red line is 
obtained when the input state is the squeezed vacuum 
state, i.e. = P(ns,0); and the black dashed line 

represents the value of the discord at the P-classicality 
threshold, i.e. = P(ng,nt). For these limiting cases, 
analytical expressions of the discord in terms of the total 
energy N of are available, even if quite cumbersome, and 
hence not reported. Being functions of a single quantity, 
they are suitable for comparison and have been plotted 
in the lower panel of Fig. 3, together with the relative 
values of the classical correlations C = Im — P- More¬ 
over, in Fig. 3 we also show, by a dashed magenta line, 
the curve corresponding to the minimum value attained 


D 



FIG. 3: (Color online) Nonclassicality by mixing with the vac¬ 
uum: In the upper panel we show the discord P as a function 
of the number of squeezed and thermal photons Us, nt for a 
balanced BS r = 1/2. The dotted red line corresponds to the 
case of a squeezed vacuum state entering the beam splitter. 
The solid blue line corresponds to a thermal input state, while 
the dashed black curve is the discord at the P-classicality 
threshold. Finally the dashed magenta curve points out the 
minimum value of the discord, obtained via numerical maxi¬ 
mization. The lower left panel shows the discord as a function 
of the total energy N-, solid blue for thermal state, dashed 
red for squeezed vacuum, black dot-dashed for P-classicality 
threshold. The right panel shows the classical correlations for 
the same input states and with the same color codes. 

by the discord (for fixed nt) obtained via numerical min¬ 
imization. 

From the left bottom panel in Fig. 3, we can see that 
the discord is a monotonically increasing function of the 
total energy N. The discord saturates to a finite value 
both for a thermal input state, for which we find 

lim P**' = log2, 

N—¥(X) 

and at the P-classicality threshold [46], where 

lim log ( 3 -^ 272 ) - ^ logy2« 0.2067. 

N—*-00 2 V / 2 

(24) 

Again we see that, for a fixed value of the total energy, 
a thermal input state results in more quantum correla¬ 
tions than a state lying on the non-separability boundary, 
although in the latter case squeezing is involved. Actu¬ 
ally, as we can see from Fig. 3, the states corresponding 
to the P-classicality threshold do not correspond to the 
states with minimum output discord (dashed magenta 
curve), conhrming again the inequivalence between P- 
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and C-classicality. The minimum output discord curve 
has been obtained numerically and we have not found any 
clear physical picture of the class of states for which this 
minimum is achieved. On the other hand, for a squeezed 
vacuum state the discord grows logarithmically for large 
N values. As a final remark, from the right plot in the 
lower panel of Fig. 3, it is apparent that as the energy 
increases the classical correlations always increase indefi¬ 
nitely, whereas as said, below the P-classicality threshold 
discord is bounded. This behavior will become clear in 
the next section, where entanglement will be considered. 

If we now release the restriction of a balanced BS and 
inquire the behavior of the discord with respect to r, we 
see that I?i| 2 (i?) and 'D 2 \i{R) differ from each others; for 
a generic value r of the transmissivity, &i-discord and 
62 -discord are simply related by an exchange of the sym- 
plectic invariants of Ii and I 2 , which amounts to a swap 
of the BS transmissivity from r to 1 — r. In Fig. 4, the 
61 -discord has been plotted as a function of r, for the 
relevant cases already mentioned. Being obtained by the 
exchange of transmissivity and reflectivity, the 62 -discord 
is simply given by a reflection about the axis t = 5 . Of 
course, in the limiting cases of transmissivity 0 and 1 , the 
discord falls to zero. Furthermore, apart for a squeezed 
vacuum input state, the behavior of I?i| 2 ( 7 ?) is not sym¬ 
metric with respect to r = ^. By increasing the incoming 
energy, the maximum of I?i| 2 ( 7 ?) ( I? 2 |i(-R) )> and hence 
the optima transmissivity, shifts towards r = 1 ( 0 ). 


Di\2 



FIG. 4: (Color online) Nonclassicality by mixing with the vac¬ 
uum: plot of the discord ©i 12 (A) as a function of the trans¬ 
missivity T for a given total energy N = 10. Colors as in Fig. 

3 . 


C. Generation of Ganssian entanglement 

Let us now focus on the generation of Gaussian en¬ 
tanglement. The explicit expressions of the symplectic 


invariants are given by 

Ii = nt{l + nt)T + N{1 - t)t , (25a) 

/2 = ^ + nt{l + nt)(l - t) -I- A^(l - t)t , (25b) 

I 3 = - [{I - nt)nt + N] {1 - t)t , (25c) 

/4 = ^ (1 + 2ntf . (25d) 


By solving the equation A_(ns, nt,T) = ^ with respect to 
rig, one finds an analytic expression for the the number 
of squeezing photons at the separability threshold: 


which does not depend on r and, most important, it 
equals the P-classicality threshold Eq. (22). This is in 
agreement with the general fact that P-nonclassicality 
is necessary and sufficient for the generation of entan¬ 
glement at a BS, regardless the Gaussian nature of the 
input state [37-44]. 

We are now going to analyze more in details the rela¬ 
tionship between the generation of Gaussian discord and 
entanglement. Since, although analytical, the expression 
of the Gaussian discord is far too involved, being in par¬ 
ticular non invertible, we proceed in our analysis by ran¬ 
domly sampling a large number of input Gaussian states 
and making them evolve trough a BS of random transmis¬ 
sivity r € [ 0 , 1 ]; for each of them, minimum symplectic 
eigenvalue of the partially transposed CM and Gaussian 
discord are then computed. 

The results are shown in Fig. 5 (light gray points), 
together with the plot obtained for evolutions trough a 
balanced BS (dark gray points). Inspecting the latter 
distribution it is easy to recover all the features already 
addressed in Fig. 3. In particular, since for input ther¬ 
mal states of large energy the discord has been found 
to reach the limiting value log 2 « 0.6931, the distribu¬ 
tion displays an asymptote, so that we can conclude that 
the region of high A corresponds to highly excited input 
thermal states. Moreover, as it can be seen following the 
dashed black line in Fig. 3, by moving on the separability 
threshold, i.e. considering the points laying on the line 
A_ = 1/2, we move from zero discord to the asymptotic 
value of Eq. (24), which is obtained for infinite input 
energy and pointed out by the black arrow. It is also 
possible to note that the minimum value of the discord 
is attained slightly below A_ = 1/2, as shown by the 
dashed magenta line of Fig. 3. 

More in general, considering arbitrary transmissivity 
(light gray points in Fig. 5) if the evolved state has a dis¬ 
cord T>ii 2 (R) > 1, it will be necessarily entangled [28, 29]: 
the avoided region of the plane {(V, A_) ] I? > 1 , A_ > 
1 / 2 } shows that the discord for separable states is al¬ 
ways smaller than one. It means that for separable— 
and hence P-classical—input states, whatever the trans¬ 
missivity, the discord between the output modes cannot 
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FIG. 5: (Color online) Nonclassicality by mixing with the 
vacuum: plot of the minimum symplectic eigenvalue of the 
partial transpose CM A_ versus the discord T>i\ 2 {R) for ran¬ 
domly generated input states Q{ns, nt) evolving through a bal¬ 
anced BS (dark gray points), and for random values of the 
transmissivity r (light gray points). The vertical dashed lines 
correspond to 1 — log 2 ~ 0.3069, log 2 ~ 0.6931, which is an 
asymptotic value for thermal states entering a balance BS, 
and 1, beyond which only entangled states (A_ < 1/2) can be 
found. The black circle stresses the portion of the plane oc¬ 
cupied by states attaining the maximum value of the discord, 
still being separable; the black arrow points the maximum 
value of the discord achieved at the separability threshold A_ 
= 1/2 in the case of a balanced BS. 


grow indefinitely [29], by simply pumping more energy. 
On the other hand, in the region 0 < I?i| 2 (i?) < 1, 
both entangled and separable states are present. An¬ 
other region of interest is the entangled region, namely 
{(P, A_) I A_ < 1/2}, where the random generated 
points get “horn-shaped”. One remarkable feature is 
that, when A_ approaches zero, i.e. entanglement is high, 
the discord becomes, loosely speaking, nearly a function 
of A_, and hence of entanglement itself. In this case, we 
note that for I? > 1 the extent of the region in Fig. 5 is 
bounded by two convergent quantities. To better clarify 
this point we set r = 1 /2 and consider an input squeezed 
vacuum state. For large iV ^ 1 the analytic expressions 
of the discord reads (at the leading order): 


superscript sq, we may write 

A_ < A)]^ « for V>1 (29) 

8 

Analogously, our numerical analysis in the same re¬ 
gion shows that the random generated points are always 
bounded from below by A™, whose expression is 

e-® 

A_ > A™ « - for V>1. (30) 

4 ^ 

Putting together Eqs. (29) and (30) we conclude that, 
for a fixed value of the discord I? > 1, the distribution of 
minimum symplectic eigenvalues of the partially trans¬ 
pose CM, is constrained in the range 

A™ < A_ < A5 . (31) 

Therefore, for I? ^ 1, A_ is an exponentially decreasing 
function of the Gaussian discord. 

Particularly interesting is finally the region corre¬ 
sponding to highly discordant—yet separable—states, 
stressed by a circle in Fig. 5. These are states sharing 
the maximum amount of quantum correlations without 
invoking entanglement. We found that these states are 
obtained for high input energies, whose value can also 
be due uniquely to thermal photons, entering a BS of 
extremely high transmissivity, namely r = 1 — e. Hav¬ 
ing unlimited thermal resources at disposal, we can still 
generate quantum correlated output states up to a value 
I?i |2 = 1, by sending a very excited input state in an 
unbalanced BS of very high transmissivity. If, always 
keeping the BS transmissivity close to one, the fraction 
of squeezed photons is such to render the input state P- 
nonclassical one, the corresponding points in the plane 
will lie just below the separability threshold, but the 
value of the discord can increase only up to the value 
1 — log2 « 0.3069 (indicated by a red dashed line). By 
further increasing the amount of squeezing, the resulting 
states will eventually occupy more entangled and more 
discordant regions of the lower branch. 


log 



+ 1 5 


(27) 


V. NONCLASSICALITY ARISING FROM 
MIXING A GAUSSIAN STATE WITH A 
THERMAL STATE 


while the minimum symplectic eigenvalue of the partial 
transpose Al*^ is 




1 

4/A 


(28) 


respectively. Form the previous equations it follows that 
in this limit Al'^ « /8, valid for high discord value. 

It turns out to be an upper bound for the random distri¬ 
bution of symplectic eigenvalues in the entangled region, 
and hence will be denoted as . Upon omitting the 


We now consider the general case, allowing for ther¬ 
mal photons to enter the second port of the BS. In fact, 
in practical scenarios a certain amount of thermal noise 
{e.g., in the form of black body radiation or scattered 
light) unavoidably participates in the interference phe¬ 
nomenon and affects the statistics of the outgoing fields. 

Though P-classicality, as already discussed, retains the 
expression (22) for the threshold value, the presence of 
another source of photons affects the properties of the 
output state. The relevant changes both to quantum 
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FIG. 6: (Color online) Nonclassicality by mixing with a ther¬ 
mal state: contour plots of Gaussian discord as a function of 
squeezed and thermal photons Us and nt in mode ai. From 
panel (a) to (d) the number thermal photons in mode a 2 are 
given by 712 = 0, 0.1,1, 3 [panel (a) is in fact the contour plot 
of Fig. 3]. Colors go from black (0.) to white (0.7). 


discord and entanglement can be again evaluated via the 
symplectic invariants, which now reads 

Ji -b 712(1 -t)'^ + T [nt + Us (1 -b 2nt) (1 - t) + ttIj] 
-b 772(1 — t) [1 -b 2ni t] , 

h + tit (1 - tY + r [712 -b Us (1 + 2712 ) (1 - r) -b tti^] 
-b 7lt(l — r) [1 -b 2 {Us + 712 -b 27ls7l2) t] , 

(32) 

h = [nl + Tit - T7s(l + 27it) - 2711712 ] (1 - t)t , 

.^4 =T^ (1 + 27lt) (1-b 2712 ) , 

16 

whence, recalling the expression for the total energy in 
mode di Eq. (3), we can see that Ii and I 2 are each 
other related via the exchange of the number of thermal 
photons Tit and 712 . 


contours of quantum discord as a function of Us and nt 
for different thermal-photon number 712 of mode 02 . We 
can see that, despite the P-classicality of the output state 
remains invariant, the C-classicality is much affected by 
the presence of an additional source of thermal photons. 
In particular, for a low number of thermal photons 712 
the region with minimal discord (darker areas in Fig. 6) 
localizes close to the P-classicality threshold, whereas it 
tends to get closer to the zero squeezed-photon axis for 
larger 712 . Also in this case, the inequivalence of the two 
notions of non-classicality is apparent. 

Since there is no threshold for the production of dis¬ 
cord, it is legitimate to enquire the generation of quantum 
correlations at a BS for the ’’cheapest” conceivable sce¬ 
nario, namely having at disposal only thermal resources 
in input. Given a certain amount of total thermal pho¬ 
tons AT = Til -b 712, which is the most convenient redistri¬ 
bution of the total energy between the two modes, in or¬ 
der to maximize the Gaussian discord at the output? The 
answer to this question is shown in In Fig. 7, where &i- 
discord as a function of the photon imbalance d = tii — 712 
between the two input modes has been plotted. For each 
transmissivity, the 61 -discord is a monotonically increas¬ 
ing function of the imbalance d, so we can conclude that 
the optimal configuration is the most asymmetric one, 
where all the thermal photons are sent in on channel, 
leaving the other in the vacuum state. Moreover, an even 
distribution of photons (d = 0 ) between the input modes 
always leads to zero output discord. This fact is appar¬ 
ent by looking at Eq. (7) where for equal input states, no 
matter the transmissivity, the correlation terms c± of the 
CM identically vanish: the phenomenon is referred to as 
transparency, since the evolution trough the BS does not 
leave any imprint on the input states. Being the optimal 
configuration the one with a thermal input in one port 
of the BS and the vacuum in the other, we already know 
that, for a given amount of energy N, there will be an 
optimal value of the transmissivity maximizing the &i- 
discord (as shown by the blue curve of Fig. 4); this fact 
is also manifest in the crossing of the blue curve and the 
red dot-dashed one in Fig. 7, corresponding to r = 0.5 
and T = 0.8 respectively, when approaching the maxi¬ 
mum imbalance. Finally, the dashed curve corresponds 
to an extremely unbalanced BS, namely r = 0.99, and, 
provided that high-enough thermal energy is available, 
the corresponding value of the 61 -discord close to the im¬ 
balance would be the greatest, eventually achieving the 
limiting value of I?i |2 = 1 , as discussed above for the 
circled points in Fig. 5. 


B. Generation of Gaussian entanglement 


A. Generation of Gaussian discord 

As in the previous Section, we first focus on the case 
of a balanced BS (r = 1/2). In Fig. 6 we show the 


As before, the equation A(7is, rij, 712 , r) = 1/2, if solved 
with respect to the squeezed number of photon Ug , gives 
a threshold on the generation of entanglement when 712 
thermal photons enter the BS. The explicit expression of 
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FIG. 7: (Color online) Nonclassicality by mixing with a ther¬ 
mal state: logarithmic plot of the 6i-discord 'D^ 2 {R) as a 
function of the imbalance d for different values of the trans¬ 
missivity r and fixed total energy N = 5. The solid blue curve 
is for a balanced BS r = 0.5, the dot-dashed red line is for 
r = 0.8 while the black dashed line is for r = 0.99. 



FIG. 8: (Color online) Nonclassicality by mixing with a ther¬ 
mal state: plot of the separability thresholds n'f'^ for fixed 
r = I and different values of n 2 . The solid blue line cor¬ 
responds to n 2 = 0, and indeed coincides with the non¬ 
classicality threshold n"^. The red dot-dashed lines repre¬ 
sents nf’’ for n 2 = 0.1, and finally the black dashed one is 
for n 2 = 1. The black lines are curves of fixed energy in 
the mode hi ni = Us -\- nt + 2nsnt : dashed for odd values 
ni = 2k + 1, k = 0,1 ,... and dot-dashed for even values 
ni = 2k, fe = 1, 2,. .. of the total energy. 


is given by 


nT = , (33) 

r(l - r) 

where Qk,i = n-kni + — (n^ — ni)T and /ri ^2 are the 

purities of the two input states. 


Contrary to the vacuum case, it is apparent that in the 
presence of a thermal state the separability threshold 
and the non-classicality threshold are no longer co¬ 
incident. In Fig. 8 several separability thresholds are 
shown, for different values of n 2 - As soon as n 2 dif¬ 
fers from zero, the number of squeezed photons re¬ 
quired to have entanglement increases — as shown both 
in Fig. 8 and Fig. 9. The previous symmetry between 

rh 

n^ 



FIG. 9: (Color online) Nonclassicality by mixing with a ther¬ 
mal state: plot of the squeezed fraction of photons Us/ni at 
the separability threshold, as a function of the total number 
of photons ni = Us + rit + 2nsnt entering the first port of the 
beam splitter, for different values of n 2 .The solid blue line 
corresponds to n 2 = 0, The red dot-dashed line corresponds 
to n 2 = 0.1, while the black dashed one to n 2 = 1. 


the notions of non-classicality in the phase space and 
non-separability no longer holds: there exists P-singular 
input states of the electromagnetic field, and hence P- 
singular output states, which nevertheless are not entan¬ 
gled. We can conclude that a hierarchy of non-classicality 
has settled down: non-separability at the output imposes 
a stricter notion of quantumness than the one put for¬ 
ward by P-singular distributions. Injecting into the BS 
a non-classical state is no longer a sufficient condition to 
get entanglement between the output modes. 

In order to better investigate this point, we express the 
separability threshold relative to n 2 thermal photons as a 
function of the P-nonclassicality threshold n^. We focus 
on the optimal case of a balanced BS, obtaining: 


sep ^ [n2 + ^«)(l-b2n2)]^ 

(l + 2n2)[l + 2h(nP)] ’ 


(34) 


where h{n^) = 71^-1- + 71^) is a monotonically in¬ 

creasing function of the non-classicality threshold. Even 
if the two thresholds and 77^ now differ, their knowl¬ 
edge enables one, given a known amount of thermal noise 
in 0 , 2 , to estimate the effective P-nonclassicality required 
in oi, i.e. how much squeezing pump into the BS, in 
order to get entanglement. When we have a squeezed 
vacuum state entering the BS in the mode oi, namely 
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rit = 0, and n 2 thermal photons in b, the separability 
threshold in Eq. (33) reduces to n|/(l + 2n2), indepen¬ 
dently on T. It is the value of the curves at n* = 0, 
as it can be seen from Fig. 8, and moreover it is the same 
expression as with n* replaced by n 2 ■ Thus, having a 
squeezed vacuum state in mode oi and a thermal state in 
mode 02 (characterized by squeezed and n 2 thermal 
photons, respectively) is equivalent to have a single mode 
Gaussian state g{ns,n 2 ) in di and the vacuum in 0 , 2 . 

Finally, in Fig. 10 we propose the same random plot 
as in Fig. 5, with the difference that random number of 
thermal photons is added in the second mode. We can 
see that the lower branch is substantially unchanged by 
the presence of thermal noise; even if the entanglement 
sets later, i.e. for higher amount of squeezing, the re¬ 
lationship with Gaussian discord remains the same. On 
the other hand, in the remaining accessible region of the 
plane, compared with Fig. 5, the points are scattered and 
the sharp pattern is now washed out. On average the dis¬ 
tribution drops towards lower values of the discord. 



0.0 0.5 1.0 1.5 2.0 2.5 


FIG. 10: (Color online) Nonclassicality by mixing with a ther¬ 
mal state: symplectic eigenvalue of the partial transpose A_ 
versus discord for randomly generated input states Q(ns,nt), 
random values of the transmissivity r and randomly gener¬ 
ated thermal states v{n 2 ) at the second port of the beam 
splitter. 


C. Effective non-classicality and non-classical depth 

As said above, in the case in which thermal photons 
are injected in the second port of the BS, the P-non- 
classicality is no longer a necessary and sufficient condi¬ 
tion to obtain output entanglement. However, remark¬ 
ably, a quantitative relation between these two notions 
can still be worked out. In particular, we will now see 
that the non-classical depth at the input determines the 
potential of generating entanglement at the output. 

Let us first consider the implicit equation defining the 
separability threshold A_(ns,n*, 712 ,r) = 1 / 2 , and let us 
call £g{T) its solution with respect to the number of ther¬ 
mal photons in the second port as a function of r. It ex¬ 


presses (as a function of the input parameters n^, nt) the 
number of thermal photons that can enter a BS of trans¬ 
missivity T in the 0,2 mode and yield an output entangled 
state. The explicitly form of ^^(t), although analytical, 
is quite cumbersome and hence has not been reported. If 
we now perform a maximization over the transmissivity 
r we obtain the following quantity: 

Eg = max Egir) . (35) 


We shall refer to this quantity as to the effective non¬ 
classicality of the state g and it embodies the maximum 
allowed number of thermal photons that can be mixed 
with g sX a. BS and still get an entangled output state. 

While the non-classical depth is a property of a single¬ 
mode state of the field, the effective non-classicality is a 
property of the two-mode configuration that we are con¬ 
sidering. In other words, the effective non-classicality Eg 
must be intended as an attempt to characterize opera¬ 
tionally the non-classical feature of a state. Thus, the 
relation between Eg and Tm, if any, is a priori unclear. 
Let us stress that the operational interpretation com¬ 
monly associated to the T^, of a single-mode state g is 
that it gives the number of thermal photons that have to 
be statistically mixed with the state g in order to obtain 
a classical state. In this sense, this operational interpre¬ 
tation of Tjn exclusively refers to single-mode states. 

In order to clarify the relation between Eg and T^ 
we notice first that, by means of a numerical maximiza¬ 
tion, it is possible to show that Eg is always obtained 
for r = 1/2. Thus the balanced BS represents the over¬ 
all optimal configuration, and in this case the effective 
non-classicality reads 




Us-nt-^ \/ns{l -I- Us) 
1 -\- 


(36) 


If now we look at the expression of the non-classical depth 
Eq. (21) and insert it in Eg, after some manipulation we 
find the following relation 




1-2W 


(37) 


Thus, the two quantities Eg and T^ which, as said, are 
defined in reference to different systems, are in fact re¬ 
lated via a simple expression. In other words, this endows 
the non-classical depth with a new operational interpre¬ 
tation: the non-classical depth of a state determines, via 
Eq. (37), the maximum number of thermal photons that 
can be mixed with it at a beam splitter without destroying 
the output entanglement. 


VI. CONCLUSIONS 

The quantum-to-classical transition for a single-mode 
bosonic system may be fully characterized by the prop¬ 
erties of its Glauber-Sudarshan P-function in the phase 
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space. On the other hand, for two-mode states, quantum¬ 
ness may be recognized either by the presence of quantum 
correlations (C-nonclassicality) or in terms of its phase 
space distribution (P-nonclassicality). In this paper we 
have addressed the generation of both types of nonclas- 
sicality by the linear mixing of a single-mode Gaussian 
state with a thermal state at a beam splitter, and have 
explored in details the relationships between the nonclas- 
sical features of the single-mode input and the P- and 
C-nonclassicality of the two-mode outputs. 

We have shown that, for mixing with vacuum, a bal¬ 
anced BS is capable of generating C-nonclassicality for 
any input state, contrary to the case of P-nonclassicality. 
In addition, the C-nonclassicality increases as the input 
P-nonclassicalitity decreases. These findings clearly con¬ 
firm in a dynamical setting the inequivalence between 
these two notions of nonclassicality that was highlighted 
in Ref. [7] in a geometrical context. We confirm this in¬ 
equivalence also for mixing with a thermal state, even if 
more complex behaviors emerge. 

In addition, we have shown that input P-classicality 
and output separability single out two thresholds which 
coincide only for the case of linear mixing with the vac¬ 
uum, whereas they are connected in a non trivial way for 


linear mixing with a thermal state. In fact, P-classicality 
at the input, as quantified by the non-classical depth, 
does determine quantitatively the potential of generat¬ 
ing output entanglement. This allows us to provide a new 
operational interpretation for the non-classical depth: it 
gives the maximum number of thermal reference photons 
that can be mixed at a beam splitter without destroying 
the output entanglement. 

By reinforcing quantitatively the inequivalence be¬ 
tween P- and C-classicality, our results paves the way 
for analyzing the dynamical relationship between differ¬ 
ent types of non-classicality in more general contexts. 
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